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Abstract 

Motivated by the analogy between successive interference cancellation and iterative belief-propagation on erasure channels, 
irregular repetition slotted ALOHA (IRSA) strategies have received a lot of attention in the design of medium access control 
protocols. The IRSA schemes have been mostly analyzed for theoretical scenarios for homogenous sources, where they are shown 
to substantially improve the system performance compared to classical slotted ALOHA protocols. In this work, we consider 
generic systems where sources in different importance classes compete for a common channel. We propose a new prioritized 
IRSA algorithm and derive the probability to correctly resolve collisions for data from each source class. We then make use 
of our theoretical analysis to formulate a new optimization problem for selecting the transmission strategies of heterogenous 
sources. We optimize both the replication probability per class and the source rate per class, in such a way that the overall system 
utility is maximized. We then propose a heuristic-based algorithm for the selection of the transmission strategy, which is built 
on intrinsic characteristics of the iterative decoding methods adopted for recovering from collisions. Experimental results validate 
the accuracy of the theoretical study and show the gain of well-chosen prioritized transmission strategies for transmission of data 
from heterogenous classes over shared wireless channels. 


Index Terms 

Random MAC strategies, slotted ALOHA, prioritized transmission schemes, successive interference cancellation, bipartite 
graphs, unequal resource allocation. 


I. Introduction 

In the era of Internet of the Things, the number of devices (sensors, machine terminal devices, portable devices, etc.) that 
are simultaneously connected to the network is expected to grow very rapidly in the near future m, m. When a massive 
number of devices share the same channel resources, there is an obvious need for networks to be opportunistically designed 
with adaptive and distributed protocols. In this context, random medium access control (MAC) protocols have received a lot 
of attention since they do not require explicit coordination between wireless network users. At the same time, when different 
classes of sources compete for a common channel, as illustrated in Fig. [T] a prioritized allocation of available resources among 
sources is necessary in order to optimize the overall network utility. The adoption of prioritized random MAC strategies in 
future networks is thus desirable, creating the need for effective optimizations of multi-sources resource allocation strategies. 

The slotted ALOHA (SA) protocol has been widely considered as one effective random MAC strategy, where users randomly 
select the time slots where they transmit information. If different users select the same time slot for transmission, a packet 
collision is experienced. While collided packets were irremediably lost in early versions of SA, recent studies have shown 
that collisions can be resolved by network diversity, multiuser detection, network coding strategies E-El, or by successive 
interference cancellation (SIC) techniques El which substantially improves the system throughput. The key concept behind 
SIC is that each user might send repetitions of the same message in different slots. If two messages from two different sources 
are sent in the same time slot (i.e., if a collision is experienced), the base station (BS) might recover the messages through 
SIC if one of the collided messages has been decoded previously. The throughput gain in applying SIC to SA schemes has 
been initially formalized in contention resolution diversity slotted ALOHA (CRDSA), where each user sends its own message 
within a MAC frame and eventually a replica in a randomly selected slot El- If a message is correctly decoded (with no 
collisions), it can be used to remove the potential interference contribution caused by the replicated message. Using similar 
concepts, authors proposed a novel random access protocol that exploits SIC in a tree algorithm in Q. A further analysis of 
such SIC-based random MAC protocol is proposed in |[8|. 

An improvement of CRDSA has been proposed in E, where the author introduced the optimized transmission technique 
for irregular repetition slotted ALOHA (IRSA) algorithm. It consists in a random SA protocol where the number of replicas 
that each user sends per frame (i.e., the replication rate) is not limited to two (or to any deterministic value) and it is rather 
randomly selected according to a pre-determined transmission probability distribution. A key connection is shown between 
the SIC in IRSA and the belief-propagation (BP) decoder of erasure codes on graphs E). This has opened the possibility of 
applying theory of rateless codes (and codes on graph in general) to IRSA schemes to optimize users’ transmission strategies 
via proper selection of their transmission probability distribution Col-IIll. For example, the IRSA has been improved with 
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Figure 1. Scenario with multiple sources (cameras) communicating to a central base station. Different levels of priority could be given to the sources, with 
high- medium- and low-priority cameras, for example, which receive different shares of the network resources. 


coded slotted ALOHA schemes m, as well as frameless slotted ALOHA protocols, in which the MAC frame size is not 
a priori selected IfTTlI . lfT2l . However, to the best of our knowledge, none of the recent works on random slotted ALOHA 
protocols with SIC have considered unequal allocation of the channel resources among users. We exactly aim at filling this 
gap and extend the IRSA framework to prioritized random MAC protocols for heterogenous wireless sources. 

We propose a prioritized IRSA algorithm, in which sources are of different classes, which transmit information to a common 
BS with a random IRSA strategy to access the channel. Within a MAC frame, each source randomly selects the replicas and the 
time slots to occupy, independently from other sources. Each prioritized source class is identified by an utility function and a 
transmission strategy. The utility function is a non decreasing function of the received rate and it is such that sources in higher 
priority classes experience a larger utility score than users in lower priority classes for a given received rate. The transmission 
strategy defined by the BS is characterized by the source rate and the transmission probability distribution (also defined as 
replication probability distribution) in IRSA, which is different in each class. The transmission probability distribution drives 
the replication rate of sources in each class, hence the performance of the system. Our objective is to find the best transmission 
strategy that maximizes the expected utility over all classes. Following the analogy between SIC and theory of codes on graph 
we analytically derive the probability to correctly resolve collisions for data in each source class along with the expected 
utility per class. Our theoretical analysis studies the performance of the iterative decoding algorithms for resolving collisions 
in unequal transmission cases. It resembles the AND-OR tree asymptotic analysis of LDPC codes over erasure channels. We 
exploit intrinsic characteristics of codes on graphs codes, i.e., the convergency of the iterative decoding method, in order 
to derive decoding probability for random transmission probability distributions. Because of the analogy between collision 
recovery schemes in IRSA and iterative message-passing algorithms on graph, we note that the UEP analysis of IRSA can 
be linked to UEP studies for irregular LDPC codes ifTTl - lfT^ . However, there is a crucial difference between the UEP MAC 
strategy considered in our work and UEP rateless coding schemes. While in the latter case the code designer controls the 
output nodes (check nodes) rather than the input nodes (message nodes), this is exactly the opposite in the IRSA case, where 
the system designer controls the input nodes (source nodes) rather than the output nodes (time slot nodes). This requires a 
different analysis of the problem compared to ini, for example. 

Our new analysis is then used to find the best unequal transmission strategy among classes in terms of both the replication 
probability and source rate per class, in such a way that the expected weighted utility is maximized. The underlying intuition 
is that more important classes should correctly receive more messages than low priority classes. This is possible either by 
sending more messages from high-priority classes (i.e., by tuning the resource allocation strategy) or by sending messages 
with a transmission rate that guarantees a lower failure probability. Results validate the accuracy of the theoretical study and 
show the gain of unequal transmission strategies for heterogenous classes. They also show that the proposed method perform 
well when compared to optimal solutions computed by simulations. 

In summary, the main contributions of this work are: 

• a theoretical study of the system performance in IRSA schemes with unequal transmission strategies, which leads to the 
asymptotic message error probability per class as well as global stability conditions; 

• a new optimization problem aimed at finding the best transmission strategy in prioritized IRSA strategies, in terms of 
both source rate and replication probability per class, for a set of heterogenous classes; 

■ a solving method based on intrinsic characteristics of the iterative decoding method adopted for recovering from collisions, 
along with proper heuristics to derive random transmission strategies for different sources. 

The reminder of this paper is organized as follows. Section HI] describes the scenario under consideration, together with key 
features of IRSA schemes. The theoretical analysis for prioritized IRSA strategies is derived in Section |III| where we also 
provide simulation results to validate the theory. The optimization problem aimed at finding the best transmission strategy for 
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Figure 2. Example of realization of IRSA strategies for two sources and four time slots per MAC frame. 


Table I 

Notation adopted in the pereormance analysis. 


Name 

Description 

Ts 

time slot over which a source can access the channel and send a replica message 

Tmac 

MAC frame duration 

N — TmacITs 

transmission slots per MAC frame 

Ck 

set of sources of class k 

Lk = \Ck\ 

number of sources assigned to class k 

Rk 

number of messages correctly received in class k 

Uk{Rk) 

utility function for class k 


prioritized sources is formulated in Sec |IV] The solving method and simulation results are also provided. Finally, we conclude 
in Section IVl 


II. Framework 

In the following, we first describe the system model considered in our work, then we detail key features of the SIC technique 
and the IRSA framework. 

A. System Model 

We consider M sources communicating to a common base station, that needs to decode the messages received from the 
sources. Sources are categorized in K priority classes and we denote by the set of active sources within the source class k. 
Let Lfe = \Ck\ be the number of sources in class k with M = Without loss of generality, we assume that the source 

classes are sorted from the most important (Ci) to the least one {Ck)- 

We assume the time axis to be discretized in MAC frames of duration Tmac and we assume that at most one source packet 
is sent per source within a MAC frame. This means that at most Lk messages can be sent per MAC frame from sources of 
class k. In our system, the sources access the channel according to the IRSA protocol m. Each MAC frame of duration Tmac 
is composed of N slots of duration Ts = Tmac/N. Each slot corresponds to a transmission interval, where one message or 
several interfering messages are sent. The traffic of the network is then computed as G = M/N. Within a MAC frame, each 
source transmits I replicas of one source message, as depicted in Eig. |2] Each replica is transmitted within one time slot and 
replicas sent from the same source are allocated to different slots, which are uniformly selected at random among the N total 
available slots. The replication rate I is selected by the source at random following a transmission probability distribution. We 
denote this distribution by {Ai k}i for sources of class k, where A; ^ is the probability that a source from the class k transmits 
I replicas within the MAC frame. 

The transmission processes are handled independently by all sources. This might lead to interference on the wireless channel. 
We assume that if a time slot is selected only by one source, the BS correctly receives the message. When multiple sources 
select the same time slot for a replica transmission, a collision is experienced. The messages interfere and the information 
transmitted over the time slot cannot be immediately recovered. However, the receiver implements successive interference 
cancellation (SIC) to partially or fully resolve collisions. This is illustrated by an example in Eig. |2] In the third time slot 
of the first MAC frame, the message sent by source 1 (mi) collides with the message of source 2 (m 2 ). This means that 
the BS receives message m 2 interfered by mi, making the messages undecodable. However, thanks to SIC techniques, the 
collision might be resolved. In particular, the BS can recover mi from the time slot 1, and once mi is revealed, m 2 can be 
“cleaned” from the interference with SIC algorithms. We assume that a perfect SIC is performed and the message is recovered 
with no errors ||9|- In a more general scenario, this interfering-cancelation procedure is iterated and may permit the recovery 
of the whole set of bursts transmitted within the same MAC frame. We refer readers to 121 for a detailed description of SIC 
techniques applied to IRSA strategies. 

Einally, we denote by i?fc G [0,..., Lk] the number of messages from class k that are correctly received at the decoder. The 
reception of Rk messages leads to an utility function Uk(Rk), which is a non-decreasing function of the rate. Let denote by 
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(a) Time slot representation (b) Graph-based representation 

Figure 3. Example of bipartite graph associated with the IRSA scheme, with 3 sources of one message each. 



Figure 4. Graph representation for the considered MAC protocol. Sources of different classes encode the information into messages that are sent in some 
of the N time slots selected at random. 


R = [i?i, i? 2 , • ■ • j Rk] the vector of received rates (or received messages) for all classes. The overall utility function for the 
system is given by 

K 

U{R) = Y,WkUk{Rk) (1) 

where Wk is a priority score that characterizes the importance of class k, and it is such that J^k ~ 1- Usually important 
classes have large weight Wk- The notation adopted in this work is summarized in Table |I] 


B. Graph-based representation of IRSA 

To study the performance of SIC strategies applied to IRSA protocols, a graph-based representation has been introduced in 
El, which shows the analogy between random MAC protocols and codes on graph. Each MAC frame can be identified by a 
bipartite graph where sources are represented by burst nodes (BNs), that correspond to messages, and transmission time slots 
are considered as slot nodes (SNs). The degree of a node defines the number of outgoing, respectively incoming edges. The 
probability of having a degree-/ BN of class k is given by A; ^ and the probability of having a degree-/ SN is denoted by ili. 
The degree of a BN corresponds to the repetition rate adopted by the corresponding source in the MAC frame. The degree 
of a SN corresponds to the number of interfering messages. Finally, the iterative message recovering procedure is associated 
to a message-passing algorithm along the graph. An example of the bipartite graph representation is shown in Fig. |2 where 
the MAC frame depicted in Fig. |3(a)| characterized by 5 time slots and 3 sources with one message each, is identified by the 
bipartite graph in Fig. |3(b)| with 5 SNs and 3 BNs. 

In Fig. m we extend the graph-based representation to the scenario with heterogenous sources that is considered in our work. 
Each source in Ck is identified by a BN class k, which transmits message replicas over the MAC frame independently from 
any other BNs. The frame status can then be represented by the bipartite graph Q = {B, S, E) where B is the set of M burst 
nodes, S is the set of N SNs, and E is the set of edges. An edge {v, e) represents a transmission of BN v in the time slot e. 
In this case, we say that u is a neighbor of e. Recall that having a degree-/ BN of class k corresponds to having a source from 
class k sending / message replicas within a MAC frame. Analogously, a degree-/ SN reflects a time slots in which / messages 
have been transmitted. The transmission strategy of the IRSA can then be identified by a node-perspective degree distribution 
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with the following polynomial representations 

Afc(x) = ^ 

i 


n{x) = ^ rtix'' 

I 


( 2 ) 


From the node-perspective degree distributions, we can also derive the edge-perspective degree distributions. We define ^ 
the probability for an edge to be incident to a degree-Z BN of class k as follows 


A;,fc = 




^^l,k 

Then, the edge-perspective degree distribution A(a;) is given by 

where AJ.(a;) = dAk{x)/dx. Analogously, the probability of having an edge attached to a degree-Z SN is 


Ifli 


Pi = 


^^i 

Hence, the edge-perspective degree distribution p{x) is given by 

p{x)=Y^pix =^. 


(3) 


(4) 


(5) 


( 6 ) 


where Vl'(x) = dft{x)/dx. 

Different {Ak{x), L^) pairs lead to different frequency of accessing the channel and different transmission rates for the 
sources of class k, creating an unequal allocation of the channels among different source classes. In the following section, we 
show how the degree distribution is used to derive a collision recovery probability per source class, which is the probability 
that a message sent from a source of a given class is successively received. 


III. Collision Resolution Probability 

We now evaluate the error probability for the prioritized IRSA schemes described above. The theoretical study is evaluated 
under the assumption of very large frame sizes (N oo), hence the analysis presented next will refer to this asymptotic setting. 
The asymptotic assumption leads to theoretical analysis which is substantially simplified but yet accurate, as already proved in 
the literature a, im. In the following, we will adopt “asymptotic setting” and “large network assumption”, interchangeably. 

In the following, first we derive the SN degree distribution in the case of prioritized transmission strategies. Then, to derive 
the decoding error probability, we extend the asymptotic analysis in a to the case of heterogenous sources. Finally, we provide 
global conditions for the stability of the iterative decoding process. 


A. Node Degree Distributions 

In our scenario, the base station assigns to each class k a transmission strategy, which is defined by a transmission distribution 
Ak{x). From Ak{x), the edge-perspective distribution Xk{x) can be evaluated as in Eq. (|4]i. The degree distribution for the 
SNs, H(x) as well as p{x), need to be computed. In a^ the SN-degree distribution is derived for the case in which all BNs 
follow the same distribution A(x). Here, we extend the analysis to the case of different distributions for different classes of 
sources, or equivalently for different types of BNs. We denote by Pk{nk) the probability that Uk edges connect Uk BNs of 
class k to the same SN. This probability is given by 





Pk = 


j^j,k 

N 


(7) 

( 8 ) 


where pk is the probability that a BN of class k has an edge incident to the considered SN. Note that pk corresponds to the 
probability that one source from class k transmits a replica in the considered time slot. The degree distribution Xli for the SNs 
then becomes 


^ Pi{ni)P2{n2)...PK{nK) 

ni,...,nK'.ni-i-...-knK=l 


s n£) 




Lk—rik 


(9) 











6 




(a) Bipartite graph. (b) Tree associated to (vi^e^,). 

Figure 5. Example of the tree representation associated to (^ 1 , 63 ) in the bipartite graph, assuming that the BN vi is of class k and it is the root of the tree 
Tj fc of depth 2i = 4. 


We can simplify Eq. (|9]l as follows. We denote by Xi the event of having the source i transmitting in a time slot t. This event 
occurs with probability pi = if source i is from the class k. Since each source independently selects the time 

slots for transmission, we have that Xi,X 2 , ■ ■ ■ ,Xm are independent Bernoulli processes, each one with its own probability 
of success Pi. Then fl; = P{Sm = 0 be modeled as a Poisson binomial process ||20l, where Sm = +X 2 + ... + Xm 

is the sum of the considered events (i.e., Sm sources transmitting in the time slot t). This permits to express Eq. as 

= X] II 

A^Fi iGA i^A'= 

where is the set of all subsets of I integers that can be selected in M, A is one of these possible subsets, and A'^ is the 
complement of A, given by {1,2,..., M} \ A. Erom 1211 . Il22l . the above expression becomes 


0/ 


M 


M 


T 


-3^ 


M 


M 


n (1 + - l)Pm) 


M- 


T 


-3^ 


j=i Lfc=i 


m—1 
■ K 


n ( 1 + {c^ - i)pk) 


Lk 


( 11 ) 


where C = exp[27rz/(M + 1)], with i = \/—l. The last equality follows from the fact that sources from the same class have 
the same probability of success pk- Einally, for large M and small pi’s, fli can be approximated by a Poisson process 

M \* / M \ / K ^ 

E Pm exp - E Pm E LkPk exp I - E LkPk 

m—1 J \ m—1 J _ \fc=l / \ k—1 


Note that we are considering large frame size networks (N -A 00 ), which leads also to large M values. Einally, each pi is 
inversely proportional to N and it becomes small for large MAC frame size (i.e., large N). This then justifies the assumptions 
considered to derive the above approximation. 

We can then derive the node degree distribution for SNs as 



M 

ri{x) = viix^ 

/=! 


exp(-x) X ^ 
/=! 


(13) 


where x = Ea^i ^kPk- Under the assumption of large networks, Eq. ( fT3] l can be further simplified as n{x) ~ exp(—x(l — x)). 
Einally, from Eq. (|6]l, the edge-perspective degree distribution for a SN becomes 


p{x) 


0'(a;) 

W) 


exp(—A(1 — x)) = exp 



AUl) 

Efc Afc 


(1-a:) 


(14) 


where G = J2k ^k/N is the traffic of the network. 
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B. Collision Recovery Probability 

We now study the collision recovery probability and analyze the asymptotic behavior of the message-passing algorithm used 
in the SIC method. We extend the analysis in ||9l and consider the case of prioritized transmission strategies. The analysis 
considers the AND-OR tree asymptotic analysis, which is frequently used in evaluating rateless codes performance ED, m 
and it has been already introduced for IRSA strategies El, in. 

First, we show the dependency among the nodes in the SIC algorithm and provide an example in Fig. |5(a)| that has two 
iterations in the decoder. We see that the slot node es receives a clean (i.e., not interfered) message from V 2 , which can then be 
decoded at the first step of the SIC procedure. We say that the edge {v 2 , 65 ) is revealed at the first iteration step and it can be 
removed from the graph, leading to a reduced graph. The decoded message from V 2 is then passed along the edge ( 1 ^ 2 , 63 ) to 
remove the interference at the third time slot of the MAC frame, namely 63 . This means that the edge {v 2 , 63 ) is also revealed 
and removed from the graph. If at the ith iteration step of the decoding the message from the burst node U 3 has also been 
decoded (and so the message along (^ 3 , 63 ) has been passed), the message from vi can be decoded through the slot node 63 . 

Dependencies between nodes in the graph can be described by a tree representation ll24ll . which permits to study the iterative 
burst decoding process. Considering a burst node v, we are interested in the probability of decoding the message through a 
slot node e. The associated tree is the one describing the neighborhood of e. This tree is rooted at v, with e as only branch 
going out from v, while the node e has branches to the neighboring burst nodes excluding v, as shown in Fig. |5] for the (ui, 63 ) 
edge. 

In more details, let denote by Ti^k the constructed tree of depth 2i with a BN of class k (e.g., ui) as a root. Each node 
at depth 2i,2i — 2,... ,2,0 are BNs (of class k for the root and of any class for the other depths value) while the nodes 
at depth 2i — 1, 2i — 3,..., 3,1 are SNs. Finally, the nodes at depth 2i,2i — 2,... ,2,0 are denoted by OR-nodes, while the 
nodes at depth 27 — 1, 2z — 3,..., 3,1 are AND-nodes in the tree representation. In the tree representation, a BN at depth 2i 
(a SN at depth 2i — 1) is marked with 1 if it is decoded (if it receives a clean message) at the ith iteration step, and marked 
with 0 otherwise. For the message sent from vi to be decoded at the ith decoding step through the slot node 63 , the slot 
node 63 has to be marked with 1. This happens if all the other {I — 1) neighboring BNs of 63 are decoded or marked with 1 
(AAD-operator), where I is the degree of the slot node under consideration. We denote by Zi the probability of the slot nodes 
at depth 27 — 1 to be marked with 0. We consider now the probability that any BN of degree I at depth i in the tree is marked 
with 1. This happens if at least one of the remaining {I — 1) neighboring SNs is marked with 1, i.e., if it has a non-interfered 
message (OR-operator). We denote by yi^k the probability that the message along the considered edge is not decoded. In our 
IRSA scheme, the messages transmitted by the sources are not known a priori and they need to be decoded. This means that 
none of the BNs is known a priori before the decoding process starts, i.e., 7/0,fc = 1,VA: and zq = 1 and all leafs are initially 
marked with 0 . 

We now evaluate the probability of having a BN of class k unknown after the 7 th iteration of the decoding process. This is 
given by 

N-l 

yi,k = F’{all AND nodes at depth 27—1 are marked with 0} = z^zl^tk = ^k{zi-i) (15) 

1=1 

where Zi-i is the probability of having a AND-node at depth 27 — 1 that is marked with 0 and A; fc is the probability for an 
edge to be incident to a degree-Z BN of class k. Each child of a AND-node is a OR-node of class k with probability qk, given 
by 

Ef=ii.A'(l) 

which is the ratio between the average number of edges going out from all BNs of class k, namely LfcA'^(l), and the average 
number of total edges in the graph, namely AfcA'^(l). It follows that 


Zi-l = 1 

= 1 


= 1 


= 1 


P{all OR nodes at depth 2? — 2 are marked with 1} 

M-1 

E 

1 — P{one OR node at depth 27 — 2 is marked with 0} 



Efc = l Qkyi—l,k 

M-l 

E 

Z=1 

K 

1 ^ ^ Qkyi—l,k 

A,-l 

1-1 

Pi 


p (1 - X] 

V k=l 


1-1 

Pi 


(16a) 

(16b) 


(16c) 

(16d) 
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Figure 6. Feedback loop associated to the iterative decoding probability described in Eq. Gl). Each eiTor probability ^ is evaluated as a function of the 
Zi—i, while Zi—i is evaluated as a function of a weighted sum of all yi—i^k probabilities. This recursion creates the feedback loop depicted in the figure. 


where in Eq. (I16bl) the probability of having one OR node at depth 2i — 2 marked with 0 is given by Qkyi-i,k- This 

corresponds to the probability of having a OR node at depth 2i — 2 and of class k marked with 0, weighted by the probability 
of having a OR node of class k as child. Substituting Eq. dlhdl i into Eq. (fTSl l. we obtain the following recursion 


yi,k — ^k {^i—1 

= 1 - p ^1 

with j/o.fc = l,Vfc, go = 1 0. Note that the error probability pi ^ is recursively derived assuming that the OR nodes at depth 
2i — 2 are actually roots of trees with depth 2i — 2 and are independent from each others. Substituting Eq. (fl4l i into Eq. (fTTl i. 
we obtain 

Einally, let Pf.{k^I) be the probability for the base station of not correctly decoding the message sent from a BN of class k 
when a iterative SIC technique is adopted with a maximum of I iterations. As already observed in Eig. |5l the message (BN) 
of class k can be decoded through any neighboring SN. It means that Pe{k, I) is computed as the probability that the message 
cannot be decoded through any edge at the tth iteration of the SIC algorithm, which means that 



) 


K 


— ^ qkyi-i,k 
k=l / 


(17) 


N 

Pe{k,I) = 


1=0 


(19) 


C. Stability Conditions 

We are now interested in evaluating the conditions under which the SIC algorithm asymptotically {i oo) converges with 
zero failure probability. We first observe that the iterative decoding process described in Eq. (fTTI i can be seen as a feedback 
loop (Eig. |6]l of which we can study the global stability, by deriving the conditions under which the system asymptotically 
converges to null error probability, i.e., to an equilibrium point (y*, z*) = (0, 0), for any initial probability {y, z). The global 
stability can be guaranteed if the error probability zt decreases at every decoding iteration, converging then to a zero error 
probability for i ^ oo. 

We first note that the control equations that characterize the feedback system, as well as the iterative decoding probability 
of Eq. (fTTl i. are 


z = l- p{i -'^Qkyk) 

k 

yk = Xk{z) 


^In the above analysis we have assume a tree ensemble representation, which implies that the bipartite graph is loop-free, since loops introduce correlation 
in the evolution of the message error probabilities. This is true for large networks. 




















Table II 

Considered polynomial distributions per class A ^. 
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Index 

Label 

Transmission Probability Distribution Ak(x) 

1 

A‘‘(x) 

0.5102x^ -1- 0.4898x"‘ 

2 

A‘’(x) 

O.SeSlx^ -1- 0.0436x^ + 0.3933x'^ 

3 

A^(x) 

0.5465a:^ -I- 0.1623a:^ -I- 0.2912a;'= 

4 

A^lx) 

a5x^Ta28aFTa22x® 

5 

A^(x) 

0.08a:^ -I- 0.14a;'‘ -|- O.Sx^ -|- O.lTx® -|- 0.14a;' -|- 0.17a:'^ 

6 

A‘(x) 

0.4977a:^ -|- 0.2207a:^ -I- O.OSSlx^ -|- 0.0756x^ -I- 0.0398a:'’+ 

0.0009a:-|- 0.0088x® -I- 0.0068x'5 + 0.0030a;i 1 + 0.0429a:i'‘ -|- 0.0081x1® ^ 0.0576x1® 



(a) BN edge-perspective 



Figure 7. Comparison of the theoretical and simulated edge-perspective degree distributions for a system with K = 2, N = 200, L\ = L 2 , and 
Ai(a;) = A^(3}) and A 2 (a^) = A^(a^). 


as shown in Fig. By substituting Eq. (fT4lt in the above equations, we obtain the following control equation for the feedback 
system 

z = l-exp(^- (20) 

'-V-^" 

f(^) 

and we can guarantee global stability if f{z) < z,\/z. This means that we should have 

1 _ exp Lc ^kLkKil) q^Xk (z)] < z . (21) 

\ k J 

In the following, we show how the global stability can be exploited to solve the transmission optimization problem in prioritized 
MAC algorithms. 

D. Analytical Performance Validation 

We now provide simulation results to validate the theoretical analysis provided above and study the performance of the 
random MAC transmission protocol in different settings. We set the maximum number of iterations used in the burst decoding 
algorithm to / = 100. For each simulated scenario, we average the experienced utility function over 1000 simulated loops. We 
test the theoretical analysis for different transmission probabilities, provided in Table |II] 

Fig.|2]first depicts the theoretical and simulated edge-perspective degree distribution for a scenario with N = 200 transmission 
slots in a MAC frame, K = 2 source classes with the same number of burst nodes (i.e., Li = L 2 ), and with the following 
degree distributions Ai(x) = A^(x) and A 2 (x) = A‘‘(a;) that have been shown to be effective in random MAC strategies ||9l- 
We see that the theoretical performance is in agreement with the simulation results, which means that the assumption of large 
networks used in the analysis still holds in the case of finite size of the MAC frame (N = 200), as observed in other works 
0 . 

We further provide a comparison between simulations and theoretical results in terms of both normalized throughput and 
utility function. We consider an illustrative scenario with two classes with priority given by wi = 0.7 and W 2 = 0.3 (this means 
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(a) Normalized Throughput for class 1 (b) Mean Utility Function 

Figure 8. Fleterogenous system performance vs. network load G for the case of K = 2, wi = 0.7, and 1 U 2 = 0.3, Li = L 2 , Ai(a;) = A‘*(a;) and 
^ 2 ( 2 :) = A“(x). Different numbers of transmission slots N are considered. 


that class 1 is more important than class 2), with the same number of burst nodes Li = L 2 , and two different distributions 
Ai(a:) = A‘^(a:) and A 2 (a;) = A^(a;). In Fig. [8] both the average of the normalized throughput and the overall utility function 
are depicted as a function of the traffic load G for 1000 simulation runs, for different size of N of the MAC frame. The traffic 
load G = ^k/N varies with the number of burst slots Li = L 2 , for a fixed N. The normalized throughput is evaluated 
as G'(l — Pe{k,I)), and we only show the throughput for class 1 since the results are mostly similar for all classes. We see 
that there is again a good match between theoretical and simulated results for all N larger than 50. When N = 50, the error 
in considering an asymptotic behavior is non negligible, but yet the utility function behavior is accurate. In particular, when 
N = 50, both theoretical and simulation values of the mean utility increases with G for G < 0.6, have a peak between [0.6, 0.8], 
and decrease beyond a traffic value of 0.8. We have obtained similar trends in other experiments for different scenarios, and 
all confirm the validity of the analysis for large MAC frames and the possibility to characterize the optimal traffic load, even 
with small MAC frames. 

We now illustrate the performance of prioritized transmission schemes for sources with different importance. In Fig. |9] the 
message error rate and the overall utility function are given as functions of the traffic load G for a sample scenario with N = 200 
transmission slots, two classes with the same number of slots {Li = L 2 ) but with different priorities {wi = 0.7 and W 2 = 0.3). 
The traffic load G = Lk/N again varies with the number of burst nodes in each class. In order to illustrate the benefits 
of prioritized transmission when sources are heterogeneous, we compare the performance of equal and unequal transmission 
probabilities for different classes. For the equal transmission strategy, denoted by EEP, we consider Ai(a::) = A 2 (x) = A®(x); 
while for the case of unequal transmission strategy, denoted by UEP, we consider different distributions Ai(a::) = A®(x) and 
h. 2 {x) = A^(x), which correspond to different transmission probabilities (namely higher replication rate for the most important 
classes). We first observe that theoretical results again match the simulations results computed over 1000 simulation runs. The 
results also show the benefit of prioritized transmission policies when source have different priorities. In Pig. |9(a)| we see 
that the message error probability per class is higher for the EEP strategy than for prioritized strategy under consideration. 
We also observe that a EEP strategy leads to a waterfall effect in proximity of a threshold value of G = 0.6, while the UEP 
strategy has a larger threshold value of G = 0.7. A larger threshold value implies a larger throughput, since a larger number 
of sources actively transmit messages within a given MAC. This gain achieved by the UEP strategy is due to the change in the 
transmission protocol for different classes: in our example, the class 2 adopts a replication rate following A 2 (a:) rather than 
Ai(a;), with a maximum node degree of 5 rather than 9. This reduces the overall number of messages sent in a MAC frame, 
hence also the load of the network. This reduced replication rate mainly affects class 2 rather than class 1, which is more 
important in our example. This leads to an overall utility function that reaches a maximum of 43.2 dB for the UEP strategy, 
as opposed to the maximum of 40.8 dB of the EEP strategy, see Pig. |9(b)| 

Another example of prioritized transmission is provided in Pig. [TO] where the same scenario of Pig. |9| is considered, but 
with a different number of burst nodes in each class, i.e., L 2 = aUi, with a = 0.1 and a = 0.33. We see again that the 
theoretical and the simulation results are generally in accordance, especially in the low and high traffic regions. We further 
see that the number of burst nodes per class affects the performance of the prioritized transmission solution. In particular, 
reducing the number of burst in the lower importance class improves the overall system performance. Pinally, we observe that 
theoretical study is an upper bound of the simulated performance and the theoretical value of G* is an upper bound of the 
simulated value of G*, where G* the threshold value beyond which the error probability rapidly reaches 1 and a waterfall 
effect is experienced. 
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(b) Mean Utility Function 


Figure 9. System performance vs. network load G for the case of N = 200, K = 2, wi = 0.7, and W 2 = 0.3, Li =1/2. For the UEP case Ai (x) = A® (x), 
and A 2 {x) = A^{x) For the EEP case, Ai(a;) = A 2 (x) = A®(a:). 



Figure 10. System performance vs. network load G for the case of UEP transmission, for N = 200, K = 2, wi = 0.7, and W 2 = 0.3, L 2 = ceLi, 
Ai(a;) = A®(x), and A 2 {x) = A*’(a;). 


Overall, the above results illustrate that prioritized transmission is beneficial when sources have different importance, and 
that the system performance is dependent on both the replication rates (or the transmission probabilities) and the number of 
burst nodes. In the following, we show how this theoretical study of MAC protocol strategies can be adopted in practical 
resource allocation optimization problems. 

IV. Prioritized Random Access Optimization 

Based on the above theoretical analysis, we formulate now a transmission policy optimization aimed at finding the best 
(Afc(x),Lfc) per class, such that the network resources are not under-utilized but not over-utilized either. More in details, 
we choose both the transmission probability and the source rate for each class, such that the overall system performance is 
maximized. The optimization, carried out by the base station, finds the best tradeoff between throughput and message decoding 
probability for the all sources. 

A. Problem Formulation 

Let denote by A = [Ai(x), A 2 (x),..., Ax(x)] the transmission policy vector with the transmission probability distributions 
for each class. By varying A and the number of bursts Lk for source class k, we can actually modify the expected received 
throughput for each class. As the throughput corresponds to an utility function score for each class, as shown in Eq. ([T]i, we 
can modify the overall system performance by tuning A and {Lfc}. In particular, we can optimize the overall utility of a system 
with heterogenous sources by maximizing the throughput of most important classes and sacrificing the throughput of least 
important ones. This leads to prioritized transmission strategies that are necessary for systems with sources that have different 
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Utility functions. Hence, instead of maximizing the overall throughput, we optimize the best transmission strategy (A, L) and 
select the one that maximizes the overall utility function, evaluated as follows 


= arg max {7(A, L) 

(AX) 


( 22 ) 


Lk 


= arg max E E -E i?2, • ■ • , RK)Pr^\Rl-,Rl, ■ • ■ , Rk',R,L) > 

{Ri=0R2=0 Rk=0 J 

where [/ (i?i, Ai,..., Rk) is the system utility function defined in Eq. ([T]i while (Pi, Pi,, Rk',R,L) is the probability 
of correctly receiving Pi, P 2 ,..., Pic messages for sources of class 1, 2,...,K, respectively. This corresponds to the 
probability of receiving Pi, P 2 ,..., Rk messages either with no collisions or messages with collisions that can be resolved 
with the iterative message recovering strategies (i.e., SIC) after a maximum number of / iterations. 

Since the transmission processes are handled independently by all sources, we can derive the probability of correctly receiving 
the the messages from different sources as follows 

K 

(23) 


pW(Pi,Pi,...,Pk;A,P) = \[pi^\Rk-Ak{x),Lk) 




where P^^\Rk; Rkix), Lk) is the probability for a source of class k to recover Rk out of Lk messages after the /-th IC 
iteration, and 1 — Pe{k, I) is the probability for the base station to correctly receive the message sent from a source of class 
k. We can then express the expected distortion as 


Li L2 Lk K Lk 

uiA,L)= E E E n 

R 1 —OR 2 —O RK—Ok—1 Rk—0 

and the problem formulation to be solved becomes 

(A*,P*) : arg max U{A,L) 

A,L 

s.t. Ak{x) < Ak+i{x) 


Lk 

Rk 


[l-Pe(fc,/)]«'=P,(fc,J) 


\Lk 


Vx G [0,1], Vfc 


(24) 

(25a) 

(25b) 


where the priority constraint in Eq. (I25bl i permits to reduce the search space. In particular, we constraint the optimization to an 
unequal recovery probability among classes such that sources from more important classes have a larger probability of correctly 
transmitting their messages compared to lower important sources. This translates in imposing Pe{k, i) < Pe{k + 1, i) for class 
k more important than class k + 1. Erom Eq. (fT9l l. the priority condition can be generalized as Akljji-i^k) < Ak+i{yi-i^k+i)- 
The optimization problem in Eq. (l25l l results in maximizing the overall system utility. 


B. Approximated Solution 

The optimization problem in Eq. ( l25l l might not be easily solved with conventional optimization frameworks. The expected 
utility function is evaluated as a weighted sum of binomial distributions, each of them having a probability Pe{k, I). Although 
Pe{k, I) can be considered as a sigmoid function to simplify the formulation, there is no a convenient optimization framework 
that is able to address the above problem jointly for both A and L variables, to the best of our knowledge. Thus, we propose a 
solving method that exploits an intrinsic property of the coded slotted ALOHA; the message error probability usually follows 
a waterfall effect 1^ . having an error probability approaching 0 for traffic network G lower than a given threshold G*, and 
rapidly approaching 1 beyond G*. The threshold value G* is usually defined as the value that is the limit of the region where 
the condition of stability hold. 

In the following, we approximate the message error probability Pe{k,I) to 0 when the stability condition is respected, i.e., 
when the convergency of the iterative message decoding algorithm is assumed. By imposing the global stability condition of 
Eq. (ITTb . we have the following instance of the optimization problem: 


K 

{A*,L*) largmax WkU{Lk) 

A.L 

s.t. Afc(a;) < Ak+i{x) 


K 


( AfcA(.(i) \ ^ 


> 1 — a; 


(26a) 

Va:G[0,l],Vfc (26b) 

Va; € [0,1]. (26c) 


fc=i 
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Algorithm 1 Prioritized Random Access Protocol Optimization 
1: step 1): Define the set £on defined as the set of pairs {A,L) in the ON region; 

£on : {(■^;^) s.t. the constraint in Eq. (I26cb is verified} 


2: step 2): Optimize the utility function within the ON region as follows 

K 

(A*,L*) : arg max > WkU(Lk) 

(A,L)g£oN 7—J 
k—1 

s.t. Afc(x) < Afc+i(x) 


Vx e [0,1],Vfc 



Figure 11. Optimized system perfonnance vs Li for the case of Ai(x) = A 2 {x) 


A®(a:). The best value of L 2 is optimized for each value Li. 


where we have assumed a null error probability if the constraint in Eq. (I26cl i is met for all classes, i.e., if the conditions for 
global stability are met. This means that the Lk messages sent from sources of class k are correctly received. The experienced 
utility function in each class is therefor Uk{Lk), and the global utility is the weighted sum of class utilities. 

In the following, we denote by ON region the set of pairs of (A,X) such that the stability condition is respected. The above 
optimization problem permits to select the ON region of the system and to seek for a solution within the region where all 
packets are decoded. This instance of our optimization problem in Eq. (l25l l offers simpler solution to the selection of the best 
transmission strategy. 

The optimization problem can be easily solved in two steps, as described in Algorithm [T] The first step (Step 1) defines the 
boundaries of the ON region by simply solving Eq. (I26cb . for the global stability. Then, finding the best pair (A*,L*) in the 
ON region reduces to solving the optimization in Step 2. This optimization has an objective function that has the form of a 
weighted sum of utility functions subject to K affine constraints. Thus, it can be easily solved for concave or linear U{Lk) 
by concave or linear programming optimization, or by more general gradient-based optimization methods for more general 
non-decreasing utility functions U{Lk). 


C. Simulation Results 

We now provide simulation results to study the performance of the prioritized random MAC transmission protocol in 
different settings. We set the maximum number of iterations used in the burst decoding algorithm to / = 100. Eor each 
simulated scenario, we average the experienced utility function over 1000 simulated loops. As utility function, we consider 
Uk{Rk) = log(i?fc),Vfc, which resembles a typical image quality metric. In Section UlI-DI we have shown that the theoretical 
study of Section IIII-BI is an upper bound of the simulated performance and the theoretical value of G* is an upper bound 
of the simulated value of G*, where G* is the end of the ON region where the error probability is negligible. While the 
asymptotic theoretical bound is surely good to design reliable replication rates A(x), it might lead to approximate solutions for 
the optimization of the resource allocation L in finite MAC size. In the following, first we show the main limitations of the 
optimization method based purely on the theoretical study. Then, we describe how these bounds combined with well-chosen 
heuristics can be adopted to jointly optimize (A, L) for selecting effective prioritized transmission strategies and we show that 
the proposed approximated algorithm still achieves performance that is close to the optimal one. 
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Figure 12. Example of ON region for a system with two source classes, N = 200, Ai(x) = A®(x), and A 2 (a:) = A“(a:). 


We first study the optimization of the number of burst nodes per class. In Fig. [TT] we show the global utility function as 
a function of the number of burst nodes in class 1 (namely, Li normalized by N), for the case of two source classes, with 
different importance wi = 0.7, W 2 = 0.3, identical transmission probability distributions Ai(a;) = A 2 (a:), and N = 100 or 
200. For each value of Li, we evaluate the best value of L 2 and plot the corresponding global utility score. The Algorithm 
[T] is considered for the optimization, where the global stability condition has been imposed by Eq. (I26cb in the first step of 
the optimization algorithm. These results are compared with simulations results where the optimal L 2 for each Li is found by 
exhaustive search. For both N = 100 and N = 200, we notice that the global stability condition imposes a stringent condition, 
leading to a tight bound of the ON region compared to the simulated one. This leads to a mean utility score that is almost 
as good as the optimal one. However, we notice that the model is not highly reliable for large values of Li. For example, 
for Li/N = 0.7 in both cases of = 100 or 200 slots, the theoretical optimization leads to a substantial drop in the utility 
function. This drop is due to the mismatch between theoretical and simulation results around G*, as discussed in Section UlI-DI 

We now better study the effect of the mismatch between theoretical stability conditions and actual ones in the optimization 
algorithm. We recall that the ON region is the one that satisfies the global stability in Eq. ( I26cb . We show in Eig. [12] the 
ON region for a case with K = 2 classes, N = 200 transmission slots, and the transmission probability distributions are 
Ai(x) = A®(a;), and h. 2 {x) = A“(a;). The ON region boundaries are derived as a function of the number of burst nodes (or 
equivalently the traffic) both from simulations and from the theoretical analysis through the global stability condition of Eq. 
(l20b . In the simulations we evaluate the ON region as the one where the decoding error probability after SIC is lower than 
10“"^. As expected, the theory gives an ON region (in blu diamonds) that is more extended than the actual one (red points), 
since the theoretical ON region gives an upper bound on the value of traffic G* that represents the transition between ON 
and OEE regions of the SIC algorithm. Unfortunately, the best L derived from Algorithm [T] approaches the boundary of the 
ON region defined by theoretical stability conditions, which is exactly the unsafe region where the theory does not necessary 
match the actual behavior of the system. This means that for most cases, the optimization of Algorithm [1] would select a 
transmission strategy such that the network is actually overloaded, which results in a poor mean utility function. Based on 
these observations, we can actually overcome the main limitation of the theoretical study with the effective following heuristics. 

The key concept is that we would like the system to work almost at the boundaries of the ON region, but not exactely 
at the boundaries. A first solution is to use the theoretical study to evaluate the boundaries, and then run simulations in a 
neighborhood of the boundaries to predict the actual ON region. This method is however not always feasible because of 
the computational complexity in simulating the considered scenario. However, for all the scenarios considered in our work, 
we empirically observe that the theoretical ON region extends beyond the actual one by about O.IG. Thus, a good heuristic 
solution, called Algorithm 2, consists in first evaluating the theoretical bound on the ON region and then translating into a “safe 
boundary” by reducing the boundaries by 10%. Einally, we can seek for the best (A, L) within the safe boundary region only, 
where the decoding error probability is actually zero. It is worth noting that the best (A, L) is selected as the best resource 
allocation L such that the utility function is maximized and (A, L) is on the actual ON region boundaries. Note that it is much 
better to work with a traffic load that is slightly lower than the optimal one G*, rather than working at a traffic load slightly 
larger than G*. Working beyond the actual ON boundaries leads to a state of error probability that is quickly approaching 
1, such that the achieved throughput can quickly fall to zero. In the following, we illustrate this statement by comparing the 
performance of Algorithm [T] and Algorithm 2. 

We propose now experiments where we compare the performance of Algorithm 1 and Algorithm 2. To evaluate the set of 
pairs (A, A) that satisfies the stability constraints in the first step of both Algorithms, Eq. (I26cl l can be solved for example by 
differential evolution ll26l as shown in El- The best polynomial distribution can also be evaluated by numerical analysis by 
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(a) Mean utility vs Li for Ai(a:) = A 2 (a:) = A‘^(a:). 



Figure 13. Optimization performance for Algorithm 1 and Algorithm 2 for a system with N = 200, and source importance given as (ii;i, W 2 ) = (0.7, 0.3). 




(a) 2 source classes, Ai(a;) = A^(a:), and A 2 (a:) = A“(a:). (b) 4 source classes, Ak(x) = 

Figure 14. Mean utility vs. the MAC frame size N, when {L^} is optimized in the case of 2 and 4 source classes. 


fixing a maximal degree for the BNs ll27l . However, for the sake of simplicity, in the following we consider the polynomial 
distributions Afc {x) derived from 0 and provided in Table HU which have been optimized to maximize the traffic threshold G*. 
Any distribution in the table is a candidate for being assigned as transmission strategy to a given source class. However, our 
optimization can be applied to any other sets of polynomial distributions that satisfies Eq. (I26cl i. We compare the performance 
of both algorithms to an exhaustive search through all {L, A) possible pairs, which leads to optimal performance in our scenario. 
For the sake of clarity, we show the results of our joint optimization over L and A as functions of one parameter at time. 

Fig. |13(a)| depicts the mean utility function as a function of the number of burst slots in the first class Li, for Ai(x) = 
A 2 (x) = A°(x), N = 200 transmission slots, and source importance {wi,W 2 ) = (0.7,0.3). For each value of Li, the number 
of messages in the second class L 2 is optimized with the different algorithms. Analogously, in Fig. |13(b)| the optimal utility 
function is depicted for different polynomial distributions Ai(a;), where the indexes on the x-axis follow the order in Table HIl 
For each distribution Ai(a:), A 2 (a:) and {Li,L 2 ) are optimized with the algorithms under comparison. For both figures, the 
optimization based on the theoretical ON region does not match the optimal results achieved by exhaustive search. This is due 
to the mismatch between theoretical and actual ON region, as described above. However, we can observe that the optimization 
based on the safe ON region in Algorithm 2 achieves a performance that approaches the optimal one. We also observe that 
by increasing the number of BNs dedicated to the most important class (i.e., Li) the global utility function increases (Fig. 
|13(a)| l. This is expected since a more prioritized transmission strategy is offered for large values of Li. If we rather look at the 
evolution of the global utility function as a function of several transmission probability distributions in Fig. |13(b)| we notice 
an almost constant behavior. This is justified by the fact that the MAC frame can be efficiently utilized and by properly tuning 
A 2 {x),Li, and L 2 for each value of Ai (x), the system achieves a large enough throughout to be in the floor region of the 
logarithmic utility function. 
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We further evaluate the performance of the optimization method proposed in Algorithm 2 in different system settings and 
provide in Fig. |14(a)| the mean utility function as a function of the number of transmission slots N, given that two classes are 
considered with Ai(x) = K^{x), and h. 2 {x) = A°‘{x). The optimization consists in finding the best pair (Ti, L 2 ) for each value 
of N, for two different pairs of weight {wi,W 2 )- The optimal performance obtained via exhaustive search is compared to the 
one obtained with Algorithm 2. In all cases, we see that the Algorithm 2 reaches mean utility scores that are almost optimal. 
As expected, we also observe that the global utility function increases with N. Finally, we run experiments in a larger system 
with K = 4: classes with importance {wi,W 2 ,W 3 ,W 4 ) = (0.6,0.2,0.1,0.1). The overall utility function is again provided as a 
function of the MAC frame size N when the number of messages is optimized with Algorithm 2 and Exhaustive search. The 
polynomial distribution is set to Ak{x) = A^{x) for all k = [1,4]. The results in Fig. |14(b)| confirm the good match of the 
performance of our heuristic-based optimization algorithm with the optimal performance. We notice again that the achieved 
overall utility function increases with the MAC frame size N, as expected. 

In conclusion, we have shown that the SIC theory can be applied to practical optimization problems, namely resource 
allocation strategies for prioritized sources. An optimal resource allocation should be evaluated either by finite-length analysis 
(which is unfortunately not available in the literature for IRSA cases) or by simulation results (not feasible in terms of 
computational complexity). Thus, we have proposed an effective heuristic solution that is practical to use and achieves 
performance approaching the optimal one. 


V. Conclusions 

We have proposed prioritized new IRSA transmission strategies for systems with sources with different levels of importance. 
We have derived a theoretical study of the system performance in IRSA schemes with heterogenous sources and analyzed 
the asymptotic message error probability per class, as well as the global stability conditions. We have then proposed a new 
optimization problem aimed at finding the best transmission strategy in prioritized IRSA, in terms of both the replication 
probability and the source rate per class. A carefully designed heuristic-based algorithm has also been developed in order 
to optimize the transmission strategy in realistic conditions. Simulation results have validated our theoretical analysis and 
demonstrated the gain of the proposed prioritized strategy. The proposed solution is practical and yet accurate, achieving 
performance close to the optimal one. This work provides the main theoretical and practical tools for a system designer to 
optimally select transmission strategies for prioritized sources communicating to a common base station in an uncoordinated 
way. 
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